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SUMMARY 

As  a warm,  well  mixed  air  mass  flows  off  a land  surface  and  over 
a cooler  sea,  the  air  is  modified  in  a layer  near  the  surface. 

Within  this  layer  humidity  decreases  while  temperature  increases 
with  height  and  a radar  duct  is  formed.  The  non-dimensional 
parameters  governing  the  growth  of  the  modified  layer  are  derived 
and  simple  forms  are  found  for  the  increase  of  layer  height  with 
fetch  and  for  the  shapes  of  humidity  and  temperature  profiles. 

From  these  relations  the  depth  and  strength  of  the  radar  duct  are 
derived  as  functions  of  the  modified  layer  depth  and  the  differences 
in  potential  temperature  and  water  vapour  pressure  between  the  over- 
land air  mass  and  air  in  thermodynamic  equilibrium  with  the  sea- 
surface. 
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INTRODUCTION 


Radar  ranges  can  be  greatly  enhanced  by  the  formation  in  the 
atmosphere  of  temperature  and  humidity  profiles  which  cause  a decrease  with 
height  over  some  interval  of  the  radio  refractive  index,  and  the  consequent 
formation  of  a radar  duct.  Such  conditions  will  result  from  an  increase 
in  temperature  and/or  a decrease  in  humidity  with  height. 

It  would  be  desirable  to  establish  the  probability  of  occurrence 
of  radar  ducting  conditions  over  the  sea  in  areas  around  Australia  from 
available  meteorological  statistics,  and  to  be  able  to  predict  such  conditions 
from  synoptic  weather  forecasts. 

To  do  either  of  these  tasks  requires  a greater  understanding  of  the 
physical  factors  involved.. 

Radar  ducts  over  the  sea  may  be  divided  into  three  main  types: 

a.  Advective  ducts  - formed,  for  example,  when  warm  air  from 
the  land  blows  over  a cooler  sea  causing  humidity  to 
decrease  and  temperature  to  increase  with  height. 

b.  Evaporative  duct  - formed  when  high  evaporation  rates 

from  the  sea  cause  sharp,  surface  based  humidity  gradients. 

c.  Elevated  ducts,  principally  at  subsidence  inversions 
where  a sinking  dry  air  mass  is  separated  from  a lower, 
molster  and  cooler  air  mass  by  a sharp  interface. 

Advective  and  evaporation  ducts  are  both  surface  based.  The 
physical  understanding  of  evaporation  ducts  when  advective  effects  are  absent 
is  well  advanced,  but  the  prevalence  of  their  occurrence  in  waters  near 
Australia  still  needs  to  be  worked  out. 

Quantitative  prediction  of  elevated  ducts  is  a major  scientific 
problem.  However  in  tropical  waters  where  the  trade  wind  inversion  appears 
to  be  a stable  feature  of  the  atmosphere  a series  of  descriptive  measurements 
could  be  useful.  An  elevated  layer  can  also  be  formed  by  advectlon  of  dry 
air  off  a land  mass  and  over  a warmer  sea  surface.  Under  these  conditions 
buoyant  convection  will  often  mix  the  air  very  well  up  to  some  level  at  which 
there  will  be  sharp  gradients  in  temperature  and  humidity.  A fuller  dis- 
cussion of  types  of  radar  duct  may  be  found  in  Lane  and  others  (1976). 
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This  note  is  concerned  with  the  formation  of  surface-based 
advective  ducts  caused  by  advection  of  a warm,  homogeneous  air  mass  over 
a cooler  sea  surface.  It  is  restricted  to  steady  state  conditions. 
Extensive  experimental  data  are  available  for  this  situation  in  reports  of 
experiments  performed  in  1944  over  the  waters  of  Massachusetts  Bay,  USA, 
(Craig,  1946,  1949;  Kerr,  1951,  Chp  3).  A summary  of  the  data  is  presented 
in  Table  1.  Although  these  are  old  data  they  are  still  the  best  available. 
Measurements  were  obtained  of  humidity  and  temperature  profiles  and  of  the 
radar  propagation  associated  with  these  profiles.  Full  details  are  avail- 
able in  the  above  references.  This  note  is  mainly  concerned  with  the 
temperature  and  humidity  profiles  and  does  not  address  the  propagation 
problem  in  much  detail. 

Advective  ducts  of  the  type  considered  here  should  be  common  off 
the  north-west  coast  of  Australia  in  daylight  hours  in  the  cooler  half  of 
the  year.  The  prevailing  winds  are  offshore  and  the  air  will  be  dry  and 
hotter  than  the  sea  for  much  of  the  day.  The  distance  out  to  sea  to  which 
the  ducts  persist  is  unknown  but  from  the  available  data  should  be  at  least 
100  km. 


2.  THEORY 

Over  land  in  reasonably  sunny  weather  the  air  is  often  well-mixed 
up  to  a height  of  the  order  of  1 km  by  buoyant  convection.  Under  such 
conditions  the  potential  temperature,  * and  potential  dew  point,  T , 
are  constant  within  the  above  height  interval  except  for  a region  of  large 
gradients  in  a thin  region  close  to  the  ground.  Under  these  conditions  the 
wind  speed,  U , and  its  direction  are  also  approximately  constant  with  height. 
When  such  a well-mixed  air  mass  moves  over  a cooler  sea-surface,  with  constant 
temperature,  Tq  , a layer  of  depth,  h , forms  in  which  the  potential 
temperature  0 , is  decreased  and  the  potential  dew  point,  T^  , increased 
relative  to  their  overland  values.  The  growth  with  fetch,  x , of  the  modi- 
fied layer  of  height  h , depends  on  the  aerodynamic  roughness,  z^  , of  the 
sea-surface,  the  wind-speed,  and  the  change  in  buoyancy  between  the  air  over 

♦ In  the  homogeneous  air  mass  the  temperature,  T , decreases  with  height,  at 
the  adiabatic  lapse  rate  of  approximately  1 C/100  m,  because  pressure  falls 
with  height.  The  potential  temperature  0 is  constant  with  height  in  a 
homogeneous  air  mass  and  T is  related  to  8 by 

0 = T + yz  . 

2 Q 

where  y = the  adiabatic  lapse  rate  - 10  C/m  . 

A similar  relation  applies  for  dew  point. 
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the  land  and  that  In  contact  with  the  sea-surface.  That  is 

h = function  ( , u,x,  z ) 

P o 

where  g = acceleration  due  to  ,ravity, 

Ap  = difference  in  density  between  air  at  the  sea 
surface  and  air  upstream, 

p = average  air  density  . 

Equation  (1)  may  be  written  in  two  alternative  forms: 


(1) 


h/x  = function  ( * , z/x) 


SA £121 

z 

Pu 


....  (2) 


>r  = function  ( * , z /x) 

z z o 

Ptl  pu 


....  (3) 


p 


It  can  be  shown  (Krauss  1972,  p 50)  that 
(0  X - T ) 3 A e 


8p 


where  A e^  = (potential  water  vapour  pressure  upstream)  - (water  vapour 
pressure  of  air  in  thermodynamic  equilibrium  at  the  sea 
surface  temperature) . 

Note  that,  if  the  upstream  air  is  relatively  hot  and  dry,  (9i-  T ) and  A e 

Ap  ° P 

are  the  opposite  sign,  so  that  if  the  former  is  small,  may  be  negative 

and  then  air  within  the  modified  layer  will  be  less  dense  than  that  aloft, 

causing  mixing  by  bouyant  convection.  This  paper  is  only  concerned  with 

3 A e A 1 s 

cases  where  Ap/p  > 0 . Usually  — — =0.1  — — , where  AT  is  the  difference 

op  I s 

between  the  upstream  potential  dew-point  temperature  and  the  sea-surface 
temperature,  and  in  most  of  the  available  experiments  AT  is  of  order  10%  of 
(6l  - T > . 


The  sea  surface  has  very  low  aerodynamic  roughness  under  most 
conditions,  and  in  the  absence  of  advective  effects  vertical  gradients  of  all 
quantities  are  small.  For  this  reason  the  influence  of  zq  will  be  assumed 
negligible  from  now  on  and  (2)  and  (3)  simplify  to 


h/x  - function 

£JL  ' 
2 

....  (A) 

^ u J 

- 2 

and  g h/u  s function  1 

2 

....  (5) 

i 

i u J 

where  g 


gAp/p . 
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In  fig  1 values  of  h obtained  from  the  profiles  of  Craig  (1946) 
and  Kerr  (1951)  are  plotted  against  x . Note  that  even  at  a fetch  of 

100  km,  h is  only  of  the  order  of  200  m . The  main  uncertaincies  in  h 
arise  from  the  difficulty  of  determining  the  exact  point  at  which  tempera- 
ture or  humidity  values  from  within  the  modified  layer  reach  the  values  in 
the  homogeneous  outer  flow.  A second  difficulty  is  that  this  point  is 
sometimes  different  for  humiditv  and  temperature.  These  two  problems  occur 
more  often  at  longer  ranges  and  suggest  that  other  processes  apart  from 
simple  modification  by  the  cool  sea-surface  may  play  a role,  especially  as 
fetch  increases.  In  addition  turbulent  fluctuations  within  the  modified 
layer,  which  will  cause  fluctuations  in  the  edge  of  the  layer,  will  have 
larger  time  and  space  scales  at  larger  fetches  and  will  thus  be  harder  to 
average  out.  It  can  be  seen  from  fig  1 that  the  scatter  between  experi- 
mental points  is  greater  than  the  uncertainty  in  the  data. 

- ,2  2 

In  fig  2,  where  g h / u is  plotted  against  g x / u , this 

scatter  is  greatly  reduced.  The  two  circled  points  had  an  overwater  fetch 
of  only  4.8  km  and  lie  above  the  rest  of  the  data  which  has  fetches  greater 
than  12  km.  There  are  less  data  in  fig  2 than  in  fig  1,  because  wind-speed 
was  always  measured  only  at  1000  ft  (305  m) , and  in  the  omitted  experiments 
an  inversion  layer  occurred  below  300  m but  above  the  top  of  the  modified 
layer.  Under  these  conditions  the  wind  speed  at  300  m will  not  be  an  appro- 
priate value  for  scaling  purposes,  because  its  relation  to  winds  below  the 
inversion  is  complex.  Figure  2 is  presented  as  a comparison  with  fig  1. 

If  the  points  omitted  in  fig  2 are  ommitted  from  fig  1 the  scatter  in  the 

latter  is  not  reduced.  However  in  fig  2 (g  /u’)  occurs  on  both  axes.  A 

„ 2 

better  presentation  is  that  of  fig  3,  where  (h/x)  is  shown  versus  (g  x/u"). 


The  straight  line  through  the  points  in  fig  2 is 


0.014 


0.55 


From  the  error  bars  in  fig  2 the  exponent  cctuLd  lie  within  t 25%  of  0.55  but 
these  limits  are  extreme  values.  It  is  more  plausible  that  0.55  Is  accurate 
within  10%  . The  straight  line  in  fig  3 is 

h/x  = 1.4b  x 10~2  ( S-4  ] (6) 

( u*  J 


This  exponent  is  also  most  likely  accurate  within  ± 10%. 


These  are  both 


empirical  relations  obtained  by  drawing  a "best-fit"  straight  line  through 
the  data  by  eye.  Calculating  a least  squares  fit  is  not  warranted  at  this 
stage. 


For  predictions  of  radar  propagation  not  only  the  depth  of  the 
modified  layer  but  the  shapes  of  temperature  and  humidity  profiles  within  that 
layer  are  required. 

In  many  turbulent  flows  mean  profiles  are  self-preserving  or  approach 

a self-preserving  form  at  a sufficiently  large  distance  downstream  where  they 

are  no  longer  influenced  by  the  upstream  boundary  conditions.  A self- 

preserving  form  is  one  in  which  the  shape  of  a profile  remains  constant  with 

fetch  if  it  is  suitably  non-dimensionalised . The  most  obvious  scaling 

parameters  for  potential  dew  point  and  temperature  profiles  are  the  height  of 

the  modified  layer,  h , the  differences,  A , between  the  upstream  potential 

dew  point  and  the  sea-surface  temperature,  and  for  temperature  profiles,  the 

difference  (0i  - T ) . 

o 

Figure  4 humidity  data  are  presented  for  x * 16  km,  as  profiles  of 

potential  dew  point  temperature  T versus  height.  They  are  presented  in  the 

normalised  form  (T  - T ) / A T versus  z/h  . There  is  a fair  amount  of 
s o s . 2 

scatter  in  these  profiles  but  no  trend  with  x or  g x/u  was  discernable 
for  x 5 16  km.  The  dashed  line  through  the  data  is 

(T  -T  ) / AT  = (z  / h)\ 

SOS 

and  this  provides  a good  fit  through  both  the  potential  dew  point  and  potential 

temperature  profiles.  All  the  data  fall  between  bounds  given  by  (z/h)1'^3  and 
1/5 

(z/h)  . The  temperature  profiles  are  presented  in  fig  5,  for  x 5 16  km, 

and  are  very  similar  to  the  dew  point  profiles. 

Profiles  of  temperature  and  dew  point  at  shorter  fetches  are  presented 
in  fig  6 . These  lie  to  the  right  of  the  profiles  presented  in  figs  4 and  5. 
The  fetch  required  for  the  shape  of  temperature  and  humidity  profiles  to  approach 
an  equilibrium  state,  if  one  exists,  will  depend  on  the  detailed  boundary 
conditions  at  the  upstream  land-sea  interface  (eg  depth  of  the  overland  super- 
adiabatic  layer,  cliffs).  These  conditions  are  only  poorly  known  in  practice. 
The  division  chosen  here  at  16  km  is  somewhat  arbitrary  and  is  based  only  on 
an  empirical  investigation  of  the  available  data. 
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To  summarise,  it  has  been  seen  that 


h/x  = function  (g  x/u‘) 


? „ 2 -0.47 

and  empirically  h/x  = 1.46  x 10  (g  x/u  ) , for  x > 5 km,  and  that 
potential  temperature  and  potential  dew-point  profiles  for  x 5 16  km  have 
the  forms 


(0  - T ) / (6j  - T ) = (z/h)'4 
o o 

(Ts  - Tq)  / ATs  = (z/h)’4  . 


Further  experiments  are  required  to  test  the  downstream  fetch  to 
which  the  above  rules  apply.  Within  100  km  of  a coastline  the  available 
evidence  suggests  they  are  reasonably  good  and  can  be  used  as  a basis  for 
deriving  refractive  index  profiles. 


3.  PROFILES  OF  REFRACTIVE  MODULUS 


The  refractive  modulus  M(z)  is  defined  by 


where 


M(z) 

N(z) 

z 

a 


(N ( z)  - 1)106  + z/a. 

refractive  index 

height,  and 

the  radius  of  the  earth. 


In  a homogeneous  air  mass  M varies  linearly  with  height,  with  a gradient  of 
approximately  0.13  m ' (Kerr,  1951)  and  the  potential  refractive  modulus 
i>  is  constant,  so  that  M = ♦ + 0.13  z . 


From  Kerr  (1951) 


* = 0 


p + b 


e 

e 


...(7) 


where  0 is  the  potential  temperature, 

e^  is  the  potential  water  bapour  pressure 

Po  is  surface  atmospheric  pressure  , 
and  c and  b are  constants  with  c = 79°K/mb  and  b = 4800°K. 
(In  inhomogeneous  air  4>  , 0 , and  e are  function  of  height  z.) 


In  an  advective  duct  it  was  found  above  that 

0 = To  + A 0 (z/h)*4  (8) 

A0  = 0i  — T and 
1 o 

T - T = A T (z/h)'4 
SO  s 


where 
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But 


(T  - T ) - [e  - e (T  )|/s  , 

s o * ‘ p wo 


and 


AT 


s ‘pi  w o 


(T  )|/s  = Ae  /s 


so  that 

where 


e **  e (T  ) + Ae  (z/h)  4 , 

n u n r\  * 


e = upstream  potential  water  vapour  pressure. 


...(9) 


e (T  ) 
w o 


water  vapour  pressure  of  air  in  equilibrium 
witli  the  sea  surface, 


and 


s = the  slope  of  the  equilibrium  water  vapour 
pressure  versus  temperature  curve  over  a 
sea-surface . 


Equation  (8)  and  (9)  may  be  rearranged  so  that 


0 ■=  8i  - Afl 


1 - (z/h) 


in  1 


....(10) 


e = e - Ae 
P Pi  P 


1 - (z/h) 4 


Put 


Then 


1 - (z/h) 


= f , and  substitute  equation  (10)  into  equation  (7). 


0 i -AOf 


( e -Ae  f • 

p + f -Ei £ — 

Ho  Oj- AO  f 


Assuming  AO  <<  0j  , this  simplifies  to 


where 


>1>  - <t>i 


# = — 

1 01 


b Ae  AO  p 

e Le 

01  01 


1 - (z/h) 


b e 


p + *- 

'o  0] 


the  upstream  potential  refractive 


| 

t. 


v 


1 

1 


. 


modulus . 
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Note  that  for  hot  dry  air  blowing  over  a cooler  sea  surface  A0  >0  , Ae  < 0 

P 

The  relative  importance  of  humidity  and  temperature  gradients  can  be  seen 

from  the  ratio  R = (b  Ae  ) / (40  p ) . 

P o 


Assuming  p = 1 000  mbar 

o 

R = 4.8  Ae  / A 0 

P 


4.8  s AT  / A0 
s 


s 'v  0(1)  mb/  C and  if  A T 'v  0(A0),  it  can  be  seen  that  the  influence 

s 

of  humidity  and  temperature  gradients  can  be  of  the  same  order  for  advective 

ducts,  but  that  the  drier  the  initial  air  the  more  important  AT  will  be. 

s 


Now  M(z)  = j'  4>j  + S 1 -(z/h)  + 0.13z  , z < h 


where 


So  that 


' 4>!  + 0 . 1 3 z 


S = — 2-  (p  A 0 - b A e ) . 


S(z/h) ” / h + 0.13  , 


z > h 


z < h 


The  height  of  the  radar  duct,  z , is  the  height  where  — = 0 , 

D a Z 

and  from  the  above  equation  this  is  at 


.(12) 


provided  z^  < h . Duct  height  is  at  z = h , if  Zp  from  equation  (12)  is 
greater  than  h , or  in  other  words,  equation  (12)  holds  for  h > (S/0.52)  . 
Note  that  from  equation  (12)  Zp  decreases  as  h (and  fetch)  increases. 

From  equation  (6) 


2 f ' 

h = 1.40  x 10"  x V5 

u 


Therefore  equation  (12)  becomes 


(S/0.52)  1.46  x 10"  x V 


2.0.16  -0.18 


9.78  S (g/u  T’1”  x 


Note  that  Zp  decreases  slowly  with  fetch.  (The  units  in  the  equations  of 

this  section  are  S.I.  apart  from  those  for  p and  water  vapour  pressures 

o 

which  are  in  mb.) 


The  forms  of  modified  refractive  index  profiles,  using  equation  (11), 
for  S = 30 jWhich  is  a typical  value,  are  shown  in  fig  7 for  a number  of  values 
of  h . The  M profile  for  homongeneous  air  and  the  locus  of  z are  shown  as 
dashed  lines. 


It  can  be  seen  here  how  the  depth  of  the  duct,  which  equals  h for 
shorter  fetches .increases  till  h = (S/0.52),  and  subsequently  decreases  with 
fetch.  Note  that  for  z^  < z < h the  gradient  of  refractive  modulus  is  reduced 
below  that  for  a homogeneous  atmosphere  so  that  the  radar  range  to  elevated 
targets  is  thereby  increased  even  for  targets  well  above  the  duct. 


A measure  of  duct  strength  can  be  obtained  by  computing  the  number 
of  modes,  m , trapped  in  the  duct  as  if  the  refractive  modulus  were  a function 
of  height  only.  This  is  given  (Kerr,  1931)  approximately  by 


rzD 


2/2  x 10 


-3 


M(z) 


M(V 


dz  + \ , 


where  X «■  radar  wave-length. 


m has  been  calculated  for  modified  layer  depths  between  20  and  200  m,  for 

S = 10,  30  and  50,  for  A=  3 cm  and  12  cm  (corresponding  to  X band  and  S band 

radars,  respectively).  The  results  are  presented  in  figs  8 (X=  3 m)  and 
9 (X  = 12  cm)  . It  can  be  seen  that  S has  more  effect  than  h on  the  number 
of  trapped  modes,  and  that  more  are  trapped  for  X = 3 cm  than  for  X = 12  cm  . 
This  means  that  although  radar  ranges  with  a smooth  sea  surface  may  be  greater 
for  3 cm  radar  the  amount  of  interference  and  consequent  fading  will  also  be 
greater  at  3 cm.  If  the  sea  surface  is  rough,  the  propagation  in  the  duct 
at  3 cm  will  sometimes  be  weaker  than  at  12  cm.  The  in-duct  attenuation  in 
dB/km  due  to  surface  scattering  n ^surface  roughness  1 ^ . From  Craig  (19A6) 

S is  typically  between  30  and  A0  so  that  the  one  or  two  trapped  modes  for 

A = 12  cm  at  these  S values,  should  give  good  propagation.  If  m = 1 exactly, 

the  horizontal  damping  is  of  the  order  of  1 dB/km  . As  m increases,  the 
damping  decreases. 


4. 


FORECASTING  OF  ADVECTIVE  DUCTS 


From  the  previous  section  it  was  seen  that  the  important  parameter 
determining  strength  of  the  radar  duct  was  S given  by 


S = o (p.AA  - bA  e ) 
01  P 


(13) 


To  determine  whether  or  not  a surface  based  advective  duct  will 
occur  at  a particular  time  one  needs  to  know  the  wind-direction  (off-shore, 
on-shore  or  parallel  to  the  coast)  and  the  sign  o f . The  actual  value 

of  | which  determines  h(x)  is  not  quite  so  important.  From  section  2 


Ap 

P 


a e 

01 


3 Ae 
J 


(14) 


Hence  the  important  quantities  for  forcecasting  advective  radar  ducting  and 
for  climatological  studies  are  the  changes  in  potential  temperature  and 
potential  water  vapour  pressure  between  the  well -mixed  upstream  air  and  that  in 
the  thermodynamic  equilibrium  at  the  sea  surface.  The  climatology  of 
advective  radar  ducting  will  be  examined  in  a future  report. 
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Figure  4.  Non-dimensional  potential  dew  point  profiles  for  x * 16  km. 

Symbols  denote  data  from  Figures  of  Kerr  (1951)  or  Craig 
(1946)  as  described  in  footnote  to  Table  1. 
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O,  K34 ; + , K35;  1,  K14;  V,  K16;  ■ ,C31. 

Dashed  line  is  (s/h)*4. 
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Figure  6.  Non-dimensional  potential  dew  point  and  temperature  profiles 
for  12.8  ? x { 16  km. 

Temperatures:  x,  C3b;  •,  K36;  0 , K38 

Dew-Point  +,  C36;  O , K36. 

Dashed  line  is  (z/h)^*. 
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Figure  7. 


Refractive  modulus  profiles  for  S = 30,  h = 50,  100,  200  m. 

Long-dashed  line  is  profile  for  a homogeneous  atmosphere 
Short-dashed  curve  is  locus  of  z^. 
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